Signatures of chaotic and non-chaotic-like behaviour in a non-linear quantum 

oscillator through photon detection. 
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The driven non-linear dufBng oscillator is a very good, and Standard, example of a quantum 
mechanical system from which classical-like orbits can be recovered from unravellings of the master 
equation. In order to generats such trajectòries in the phase space of this oscillator in this paper 
we use a the quantum jumps unravelling together with a suitable application of the correspondence 
principle. We analyse the measured readout by considering the power spectra of photon counts 
produced by the quantum jumps. Here we show that localisation of the wave packet from the 
measurement of the oscillator by the photon detector produces a concomitant structure in the 
power spectra of the measured output. Furthermore, we demonstrate that this spectral analysis can 
be used to distinguish between different modes of the underlying dynamics of the oscillator. 



I. INTRODUCTION 



There is currently an intense interest being shown in 
the possible application of quantum devices to fields such 
as computing and Information processing Q]. The goal 
is to construct machinery which operates manifestly at 
the quantum level. In any successful development of such 
technology the role of measurement in quantum systems 
will be of central, indeed crucial, importance (see for 
example 0). In order to extend our understanding of 
this problem we have recently investigated the coupling 
together of quantum systems that, to a good approxi- 
mation, appear classical (via the correspondence limit) 
but whose underlying behaviour is strictly quantum me- 
chanical 01 ■ In this work we foUowed the cvolution of 
two coupled, and identical. quantised Dufhng oscillators 
as our example system. We utilised two unravellings of 
the master equation to describe this system: quantum 
State diffusion and quantum jumps which correspond, 
respectively, to unit-efRciency heterodyne measurement 
(or ambi-quadrature homodyne detection) and photon 
detection We demonstrated that the entanglement 
that cxists between the two oscillators depends on the 
nature of their dynamics. Explicitly, we showed that 
whilst the dynamics was chaotic-like the entanglement 
between the oscillators remained high; conversely, if the 
two oscillators entrained into a periòdic orbit the degree 
of entanglement became very small. 

With this background we subsequently became inter- 
ested in acquiring a detailed understanding of experimen- 
tal readouts of quantum chaotic-like systems. In this 
paper we have chosen to explore the subject through 
the quantum jumps unravelling of the master equa- 
tion Here, the measured output is easily iden- 
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tified, namely a click or no click in the photon detector. 
However, this measurement process is rather unique in 
the fact that it possesses no classical analogue. Indeed, 
this is the case even when the system under considcration 
may appear to be evolving along a classical trajectory. 
Interestingly, despite the fact that the photon detector 
has no classical analogue, it is the very presence of this as 
a source of decoherence that is responsible for recovering 
classical-like orbits in the ((g) , (p)) phase plane (despite 
the fact that we measure neither q nor p). The subject of 
recovering such chaotic-like dynamics from unravellings 
of the master equation has becn studicd in depth in the 
literature 0, 0, 0, lïol ITÏI l and a detailed discussion is 
beyond the scope of this paper. However, we notc that 
recently in [l^ resonances have been observed in a model 
of a non-linear nano-mechanical resonator that is absent 
in the corresponding classical model. In this present work 
we have chosen to scale the oscillator so that we recover 
orbits similar to those generated from a classical analysis. 



II. BACKGROUND 

In this work we study the output resulting from the 
measurement of quantum objects where the measurement 
device generates decoherence effects. In this limit the sys- 
tem exhibits dynamical behaviour in terms of its expec- 
tation vàlues very much like those observed in its classical 
counterpart. In this work we investigate the region of pa- 
rametcr space under which the classical system exhibits 
chaotic motion. Of the many models that could be used 
we have chosen the Quantum Jumps approach 0, S IBl ■ 
We note that this is only one of several possible unrav- 
ellings of the master equation that correspond to the 
continuous measurement of the quantum object consid- 
ered. Our motivation for using this approach is that the 
recorded output of the measurement is completely trans- 
parent i.e. the photon counter either registers a photon 
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or it does not. opcn quantum system is given by thc stochastic Itò in- 

In thc quantum jumps unravelling of the master equa- cremeiit equation 
tion thc cvolution of the (pure) state vector for an 



--H[iP)dt 



1^) dN, 



(1) 



whcrc H is the Hamihonian. Li are thc Linblad oper- 
ators that represent couphng to the environmcntal de- 
grees of freedom, dt is the time increment, and dNj is 
a Poissonian noise process such that dNjdNk = SjkdNj, 

dNjdt = and dW] = (^^Lj^j dt. These latter condi- 
tions imply that jumps occur randomly at a rate that is 
determincd by i^L^^Lj^. Wc will find that this is vcry 

important when explaining the results presented later 
in this paper. For an exceUent discussion of quantum 
trajectòries interpreted reahstic model of a sys- 

tem that is being continuously monitored see Q. For 
an interesting and more general discussion on the emer- 
gence of classical-like behaviour from quantum systems 
see (Ulli. 

The Hamiltonian for our, Standard, example system of 
thc Duffing oscillator is givcn by 

where q and p arc thc canonically conjugatc position and 
momentum operators for thc oscillator. In this example 
we have only one Linblad operator which is L = '\/2ra, 
whcrc a is thc oscillator annihilation (lowcring) operator, 
g is thc drive amplitudc and F = 0.125 quantifics thc 
damping. 

In order to apply the corrcspondcncc principal to this 
system, and recover classical-likc dynamics, we have in- 
troduced in Eq. Q the paramcter (3. For this Hamil- 
tonian it has two interpretations that are mathemati- 
cally equivalent. Firstly, it can be considered to scale 
h itsclf, or, altcrnativcly we can simply vicw (3 as scal- 
ing thc Hamiltonian, Icaving h fixcd, so that the relativc 
motion of thc cxpectation vàlues of the observables be- 
comes largc compared with the minimum area (/i/2) in 
the phase space. In either casc, the system behaves more 
classically as (3 tends to zero from its maximum value of 
one. In this work we have chosen to set j3 = 0.1. 

III. RESULTS 

Lct us now consider thc spccific example of a DufRng 
oscillator with a drive amplitude g = 0.3. This param- 



cter, together with all thc those already specified, form 
thc clàssic example used to demonstrate that chaotic-like 
behaviour can be recovered for open quantum systems by 
using unravcUings of thc master equation 0, 0,M ^3 ■ In 
Fig. ^we compare thc powcr spcctra of the classical po- 
sition coordinate with that of (g). Here noise has been 
added to the classical system so as to mimic the level of 
quantum noise that is present in the stochastic elements 
of our chosen unravelling of thc master equation and we 
have solved for a realisation of the Lagivan equation. As 
can be secn, for this value of f3 there is a very good match 
bctween these two results. Morcover, both display power 
spcctra that arc typical for oscillators in chaotic orbits. 

However, it is not position that is thc measurcd out- 
put in this model, but thc quantum jumps recorded, as a 
function JV{t) of time in thc photon detector. As stated 
above, these jumps occur randomly at a rate that is de- 
termincd by which, for this example, is 2F(n). 

Hence, the probability of making a jump is proportional 
to the number of photons in the state of the system at 
any one time. 

We now consider a spccial casc that occurs frcquently 
in thc classical limit, namcly whcrc \ip) localiscs approx- 
imately to a coherent (Gaussian) state. It is apparent 
that for such a state thc chance of observing a jump is 
proportional to the square of the distance in {{q) , (p)) of 
the state from the origin. In order to illustrate thc impli- 
cations of this, lct us consider a drivcn simple harmònic 
oscillator. Thc Hamiltonian is 

Hs = Íp2 + Íg2 + |cos(í)g 

and we note that in this spccial case the only effect of 
/? = 0.1 is to scale the amplitudc of the drive (again we 
set g = 0.3). We now solvc Eq. |^ using this Hamiltonian 
and allow thc system to scttlc into a stcady state. Thcn, 
as thc phase portrait for this system simply dcscribcs a 
circle centred about (0, 0) we would expect the power 
spcctra of photons counted to be the same as those for 
white noise. Indeed, this is clearly seen in Fig. |21 where 
we show the power spectrum for both (a) the position 
operator and (b) the measurcd quantum jumps. 

For more complicatcd orbits, such as those exhibited 
by thc Duffing oscillator, we would expect to see some 
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(a) Classical Duffing oscillator. 



(a) Power spcctrum of (q). 
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(b) Quantum Duffing oscillator. 



(b) Power spectrum olAÍ{t). 



FIG. 1: Power spectrum of the position, x for the classical 
Duffing oscillator and (q) for the quantum Duffing oscillator 
f} — 0.1. The frequency is normalised to the drive frequency 
of the oscillator. 



FIG. 2: Power spectrum of the position (g) and photons 
counted J^{t) for the quantum simple harmònic oscillator in 
a steady state. Here the frequency is normalised to the drive 
frequency of the oscillator. 



evidence of the underlying dynamical behaviour. Hcncc, 
localisation of \ip) from the measurement of the Dufhng 
oscillator through the photon detector forms a concomi- 
tant structure in the power spectrum of the mea sured 
output. In Fig. 13 we show, for comparison with l(b) 
such a power spectrum. 

As we can see from Fig.|31the power spectrum for this 
chaotic mode of operation reveals some structure. How- 
ever, it is not clear from this picture alone how wc might 
relate this result to that shown in Fig. l(b) It is therefore 
reasonable to ask if this result does indeed teli us any- 
thing about the underlying dynamics of the oscillator. 
We have addressed this point by computing the power 
spectrum of both (q) and Af(t) for drive amplitudes in 
the range < g < 3, the results of which are presented 
in Fig. 21 

Although the functional form of these power spectra 
obviously differ, they do clearly exhibit changes in be- 
haviour that are coincident in the drive amplitudes of 



both figures. These are identified as intervals in g la- 
belled I, II, . . . in Fig. H 

To help clarify Fig. ^ we provide in Fig. [S] explícit 
power spectra of both (q) and the quantum jumps J\í(t) 
for regions I-IV. As expected for region I in Figs. |5(a)| 
and 5(b) we see a strong resonance at the frequency of 
the drive. The broadband behaviour characteristic of 
the chaotic phenomena associated with region II is evi- 
dent in Fig. 
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and a concomitant, although different, 
structure in tE^ power spectrum of the detected photons. 
In region III of Fig. 0] we again return to a periòdic or- 
bit. In Fig. |5(e)| the power spectra for (q) exhibits a 
peak at the drive frequency, however, the power spectra 
of Aí{t) peaks at twice this frequency. The lack of coin- 
cidence between these two figures will be ex plai ned fu Uy 
in the foUowing text. Finally in Figs. 5(g) and 5(h) we 
see power the power spectrum of the quasi periòdic dy- 
namics of region IV, again the discrepancy between these 
two figures is discussed bclow. 
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FIG. 3: Power spectrum of the me asured quantum jumps 
jV{t) for the Duffing oscillator of Fig. [ï(b)1 



s* 

<; 
I 

c I 
Q 



~i — ^ r 



Normalisecl Frectuency 



I 

H 

Q 



(b) 



Normalised Frequency 

FIG. 4: (color online) Power spectrum of the (a) (q) and (b) 
measured quantum jumps as a function of drive amplitude. 



The mechanism through which the detection of pho- 
tons can yicld significant information about the undcrly- 
ing dynamics of the systcm can bc understood by looking 
at the phase portraits of (q) and (p) associated with the 
regions I-IV of Fig. 01 for thosc vàlues of drive used in 
Fig. |S1 these are shown in Fig. |ül 

For region I there is a strictly periòdic response on 
both power spectra at the drive frequency of the oscilla- 
tor. It can be seen from Fig. 6(a) that, because of the 
distance from the origin, the chance of there being a pho- 
ton counted at point A is morc likcly than at point B. 



As this occurs at the same frequency as the oscillations 
of ((?), we have direct agrccment in the position of the 
resonance in each of the different spectra. 



In region II, and as is clear from Fig. |6(b) the system 
is foUowing a chaotic-likc trajectory. Although the power 
spectra differ drastically in their structure, they do both 
exhibit broad band behaviour that is characteristic of 
chaotic or bits. 

As the drive amplitude is increased further, region III 
in Fig. 01 is accessed as the behaviour observed in region 
II ceases. For this range of drive amplitudes the solution 



is again a stable periòdic orbit as displayed in Fig. 6(c) 
However, this time, whilst the power spectrum of (g) ex- 
hibits a resonance at the drive frequency, that of Aí{t) 
appears at double this frequency. The explanation for 
this is simply that the probability of detccting a photon 
when the orbit is in a region of phas e spac e near the ori- 
gin, such as those marked C in Fig. |6(c)| is less that at 
those further away as in the region of D. This variation in 
probability occurs twicc a period and thereforc produces 
a resonance at double the drive frequency. An immedi- 
ate coroUary is that, by detecting a resonance at either of 
these different freqüències in the power spectra of JV{t) , 
we can determine whether the oscillator is in region I 
or III of Fig. 01 From our analysis in Q it may, in some 
circumstances, be advantageous to place the system in a 
chaotic orbit. It is possible that this sort of analysis could 
be used to increase or decrease drive amplitude as part of 
a feedback and control clement for quantum machinery. 

Finally, the power spectrum of (q) in region IV of 
Fig. 01 (a) is characteristic of quasi-periodic behaviour. 
Using a similar argument to the one above, we can trans- 
fer these features onto the spectrum of Aíit) . If we com- 
pare this result with the, albeit noisy, phase portrait 
of Fig. 6 (d) there is clcar evidence of quasi-periodic be- 
haviour. 

We have demonstrated, using the Duffing oscillator as 
our example system, that the different features exhibited 
in the power spectrum of the photon count can be asso- 
ciated with concomitant features in the power spectrum 
of the position operator (and vice versa). We note that 
for any given experimental system where there is a di- 
rect correspondence between the power spectrum of Af 
and (x) that the power spectrum of providcs us with 
the same amount of information about the undcrlying 
dynamics (e.g. chaotic. quasi-periodic ctc.) as the power 
spectrum of (x) . We would like to emphasise that if this 
direct correspondence did not exist then we would not 
necessarily be able to make such an assertion. For ex- 
ample, this situation might occur for a system in which 
there was a high dcgrce of symmetry in the (x) — (p) phase 
portrait. However, such a detailed study is bcyond the 
scope of this paper. 



5 




Normalised Prequency Normalised Frequency 

(a) Power spcctrum of (q), g = 0.1. (b) Power spectrum of Af{t), g = 0.1. 
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(c) Power spectrum of (q), g = 0.3. (d) Power spectrum ol J\í{t), g = 0.3. 
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(e) Power spectrum of (ç), g = 1.25. (f) Power spectrum of Aí{t), g = 1.25. 
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(g) Power spectrum of (q), g = 2.5. 



(h) Power spectrum of Af{t), g = 2.5. 
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FIG. 6: Example phase portraits for four difFerent drive amplitudes corresponding to the regions I to IV as marked in the 
power spectrum of Fig. 0] 
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IV. CONCLUSION 

In this work we havc shown that, via analysis of the 
power spectra of thc photons detected in a quantum 
jumps model of a Duffing oscillator, we can obtain signa- 
tures of the underlying dynamics of the oscillator. Again, 
we note that the decoherence associated with actually 
measuring these jumps is that which, through localisa- 
tion of the state vector, enables these classical-like orbits 
to become manifest. We have also demonstrated that 
the power spectra of the counted photons can be used to 
distinguish bctwccn diffcrcnt modes of operation of the 



oscillator. Hence, this or some form of time-frcquency 
analysis, could be used in the fccdback and control of 
open quantum systems, a tòpic likcly to bc of interest in 
some of the emerging quantum technologies. 
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